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Abstract. Primordial black holes (pbhs) have motivated many studies since it 
was shown that they should evaporate and produce all kinds of particles (Hawking 


1974). Recent experimental measurements of cosmic rays with great accuracy, the¬ 


oretical investigations on the possible formation mechanisms and detailed evap¬ 
oration processes have revived the interest in such astrophysical objects. This 
article aims at using the latest developments on antiproton propagation models 
(Maurin et al. 


2001 


and Donato et al. 20016) together with new data from BESS, 
CAPRICE and AMS experiments to constrain the local amount of pbh dark mat¬ 
ter. Depending on the diffusion halo parameters and on the details of emission 
mechanisms, we derive an average upper limit of the order of ~ 1.7 • 10“®®g 
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Introduction 

Very small black holes could have formed in the early universe from initial density in¬ 
homogeneities (Hawking 1971), from phase transition (Hawking 1982|), from collapse of 


cosmic strings (Hawking 1989 ) or as a result of a softening of the equation of state 
(Canuto 197^). It was also shown by Choptuik (1993) and, more recently, studied in 


the framework of double inflation (Kim 200C), that pbhs could even have formed by 
near-critical collapse in the expanding universe. 

The interest in primordial black holes was clearly revived in the last years for several 


reasons. First, new experimental data on gamma-rays (Connaughton 1998) and cosmic 
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rays (Maki et al. 1996 ) together with the construction of neutrinos detectors (Bugaev & 
Konishchev 2001), of extremely high-energy particle observatories (Barrau 2000|) and of 


gravitational waves interferometers (Nakamura et al. 1997) give interesting investigation 
means to look for indirect signatures of pbhs. Then, primordial black holes have been used 
to derive a quite stringent limit on the spectral index of the initial scalar perturbation 
power spectrum due to a dust-like phase of the cosmological expansion (Kotok & Naselsky 


1998). It was also found that pbhs are a great probe of the early universe with a varying 


gravitational constant (Carr 2001 ). Finally, significant progress have been made in the 
understanding of the evaporation mechanism itself, both at usual energies (Parikh & 


Wilczek 2000) and in the near-planckian tail of the spectrum (Barrau & Alexeyev 2001, 


Alexeyev et al. 2001 ). 

Among other cosmic rays, antiprotons are especially interesting as their secondary 
flux is both rather small (the p/p ratio is lower than 10“^ at all energies) and quite well 
known. This study aims at deriving new estimations of the pbh local density, taking into 
account the latest measurements of the antiproton spectrum and realistic models for cos¬ 
mic ray propagation. Contrarily to what has been thought in early studies on antiproton 


spectra (see, e.g., Gaisser & Schaefer 1992 for a review) and to what is still assumed in 


some recent pbh studies (e.g. Kanazawa et al. 200C ) there is no need for any exotic as- 
trophysical source to account for the measured p flux. On the other hand, the very good 
agreement between experimental data and theoretical predictions with a purely secondary 


origin of antiprotons (see Donato et al. 2001a) allows to derive quite stringent upper limit 
on the amount of pbh dark matter. This article focuses on this point and is organised 
in the following way : the first section reminds general consideration on the evaporation 
process and on the subsequent fragmentation phenomena. The second part is devoted to 
the primary antiproton sources while the third part describes all propagation aspects, 
with a particular emphasis on the uncertainties coming from astrophysical parameters 
(mostly the magnetic halo thickness) and nuclear physics cross-sections. Next, the re¬ 
sulting “top of atmosphere” (TOA) antiproton spectra for different formation conditions 
and for different possible emission models are studied. Finally, we derive the resulting 
upper limits on the local pbh density and consider possible future developments both on 
the experimental and theoretical sides. 


1. Antiproton emission from pbhs 
1.1. Hawking process 

The Hawking black hole evaporation process can be intuitively understood as a quantum 


creation of particles from the vacuum by an external field (see Frolov & Novikov 1998 
for more details). This can occur as a result of gravitation in a region where the Killing 
vector is spacelike, i.e. lying inside the = 0 surface, which is the event horizon in 

















Barrau et al.: Antiprotons from primordial black holes 


3 


a static spacetime. This basic argument shows that particle creation by a gravitational 
field, in a stationary case, is possible only if it contains a black hole. Although very 
similar to the effect of particle creation by an electric field, the Hawking process has a 
fundamental difference: since the states of negative energy are confined inside the hole, 
only one of the created particles can appear outside and reach infinity. This means that 
the classical observer has access to only a part of the total quantum system. 

To derive the accurate emission process, which mimics a Planck law. Hawking used 
the usual quantum mechanical wave equation for a collapsing object with a postcollapse 
classical curved metric (Hawking |l975|) . He found that the emission spectrum for particles 
of energy Q per unit of time t is, for each degree of freedom: 


d"iV ^_ Ts _ 

dQdt ('exp - (-1)2«) 

where contributions of electric potential and angular velocity have been neglected since 
the black hole discharges and finishes its rotation much faster than it evaporates (Gibbons 


1975 and Page 1977). k is the surface gravity, s is the spin of the emitted species and Ps 


is the absorption probability. If we introduce the Hawking temperature (one of the rare 
physics formula using all the fundamental constants) defined by 
hc^ IQiSg 


T = 


IQ'KkGM 


M 


-GeV 


the argument of the exponent becomes simply a function of Q/kT. Although the absorp¬ 
tion probability is often approximated by its relativistic limit, we took into account in 


this work its real expression for non-relativistic particles: 
iTras{Q,M,fi) 2 


P., = 


h'^c^ 




where as is the absorption cross-section computed numerically (Page [l976 ) and /r is the 
rest mass of the emitted particle. 


1.2. Hadronization 


As it was shown by MacGibbon and Webber ( |1990 ), when the black hole temperature 
is greater than the quantum chromodynamics confinement scale Aqcd^ quarks and glu¬ 
ons jets are emitted instead of composite hadrons. To evaluate the number of emitted 
antiprotons, one therefore needs to perform the following convolution: 


d^iVp 

dEdt 


= E 


/ o 


Q=E 




1 X dg,p(g,ig) ^Q 

drj 


where aj is the number of degrees of freedom, E is the antiproton energy and 
dgjp{Q,E)/dE is the normalized differential fragmentation function, i.e. the number 
of antiprotons between E and E + dE created by a parton jet of type j and energy Q 
(including decay products). The fragmentation functions have been evaluated with the 


high-energy physics frequently used event generator pythia/jetset (Tjostrand 1994 ), 
based on the so-called string fragmentation model. 
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2. Primary sources 


In order to compute the antiproton spectrum for a given local pbh density, the number 
qp^^{r, z, E) of antiprotons emitted with kinetic energy between E and E-\-dE per unit 
volume and time, must be evaluated. It is proportional to the number d^n/dMdV, of 
PBHs per unit of mass and volume and to the individual flux cPNp/dE dt emitted by one 
PBH, so that 


qP^^^{r,z,E) 


d^NpjM, E) d^n{r,z) 
dEdt ' dMdV 


where r and z are the cylindrical coordinates describing position in our Galaxy. As the 
physics of evaporation and the mass spectrum of pbh do not depend on their numerical 
density, the primary source term can be split into a spatial and a spectral dependance, 
as 


A) = gP"™(r,2) x Q^^^{E) 


with 

= and Q^^^{E) 

Pq 



d^Np 
dE dt 


d?nQ 

dMdV 


dM 


( 1 ) 

( 2 ) 


2.1. Spatial distribution 


Primordial black holes should have followed the dark matter particles during the forma¬ 
tion of halos, so that they may have the same spatial distribution. Unfortunately, this 
distribution is not very well known, and several independent evidences are contradictory. 


In the absence of a clear answer to this problem {e.g. see Ghez. et al. 1998, Gondolo & 


Silk 1999 , Debattista and Sellwood 1998 , 2000 ), several profiles for the pbh distribution 
can be used, with the generic form 


„PBH( 


",z) 


PBH 

P<3 


Rp> 




TDOi I TDCX. 
Xtp “T 


■0 


( 3 ) 


-I- + z'^) 

where spherical symmetry has been assumed. Different cases have been considered, with 
numerical values taken in Galcaneo-Roldan and Moore ( ^001 ). Numerical simulations 
point toward singular profiles with 7 = 1.5, a = 1.5, e = 1 and Rc = 33.2 kpc (Moore 
1999| ) or 7 = 1, a = 1, e = 2 and Rc = 27.7 kpc (Navarro, Frank & White 1996t hereafter 


NFW). We also considered an isothermal profile with 7 = 0, a = 2 and e = 1 and a 
modified isothermal profile with 7 = 0, a = 2 and e = 1.5 and Rc = 24.3 kpc. 


2.2. Spectral dependence 

It can be seen from eq. (H) that the keypoints to derive the Q^^^{E) term are the shape 
of the mass spectrum today and the boundaries of the integral. The choice of Mmin and 
Mmax is discussed later. 
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2.2.1. PBH mass spectrum 

The mass spectrum today is the result of the evolution of the initial mass spectrum in 
time. The usual picture is based on the idea that a PBH can form if an overdense region 


collapses a contrast density S (Harrison 197C , Carr & Hawking 1974 ) so that 1/3 < d < 1. 
For a scale invariant power spectrum, this leads to 


dn 


oc M, 


- 5/2 


This mass spectrum will be regarded as the standard model (see, e.g. Bugaev & Konishev 


2001 for a review). To deduce the mass spectrum today from this initial one, the mass loss 
rate must be evaluated: dM/dt « {7.8ds=i/2 + 3.1ds=i} x 10^^ g^ s“^ with ds=i /2 = 90 
and ds=i = 27 (MacGibbon 199l| ). If we approximate the previous law with a « const, 
which is correct as far as the number of degrees of freedom does not increase dramatically 
fast, the mass can be written as a function of time as: 

M,„,t«(3at + M3)i/3. 

It is then straightforward to see that with 
dn _ dn dM^nit 
dM dM^nit dM 

the resulting spectrum today is characterised by 
dn 
dM 
dn 
dM 

where M* « 5 x lO^'^g is the initial mass of a pbh expiring nowadays. In this model, the 
mass of a pbh formed at time t is determined by the horizon mass at this epoch, 

I Mpi 

M = - —t. 

8 tpi 

In the usual cosmological picture, only pbhs formed after inflation should be taken into 
account as those produced before would be extremely diluted due to the huge increase 
of the cosmic scale factor. The end-time of inflation tpp being related to the reheating 
temperature Tpp by 


(X M^ for M < Mh, 
oc M-5/2 for M > M* 


tpH ~ 0.35 


Mpi 

T‘2 

^RH 


where g « 100 is the number of degrees of freedom in the early universe, the minimal 
mass of PBH in the initial spectrum can be given as a function Tpp. 

It should be emphasized that the shape of the spectrum below M* « (3at)^/^ does 
not depend on any assumption about the initial mass spectrum. Whatever the initial 
spectrum, it should increase as M^ today: this is only due to the dMinu /dM term which 
is proportional to M^ for small masses whereas the dn/dMinn term is nearly constant. 
As shown in the next section, the antiproton emission is governed by pbhs with masses 
below lO^'^g. The results derived are therefore independent of the details of the formation 


mechanism. 
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Fig. 1. Primary antiproton flux with the standard mass spectrum before propagation (in 
arbitrary units). Curve (1) is for M e [Mpi, lO^^g], curve (2) is for M e [lO^^g, lO^^g], 
curve (3) is for M £ [lO^^g, 5 • lO^^g], curve (4) is for M > 5 • lO^^g and the thick line is 
the full spectrum. 

2.2.2. Cumulative source before propagation 

Fig. I gives the antiproton flux after convolution with the pbh mass spectrum before 
propagation. As it could be expected, pbhs heavier than M* nearly do not contribute as 
both their number density and their temperature is low. On the other hand, only pbhs 
with very small masses, i.e. with very high temperatures, contribute to the high energy 
tail of the antiproton spectrum. It should be noticed that the only reason why they do 
not also dominate the low energy part is that the mass spectrum behaves like below 
M*. In a Friedman universe without inflation, this important feature would make the 
accurate choice of the mass spectrum lower bound Mmin irrelevant (as long as it remains 
much smaller than M*), due to the very small number density of black holes in this mass 


region. 
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Fig. 0 shows how inflation modifies the standard mass spectrum. The reheating tem¬ 
perature used herafter is nothing else than a way of taking into account a cutoff in the 
mass spectrum due to the rather large horizon size after inflation. There is nearly no 
constraint, neither on the theoretical side nor on the observational one, on this reheating 
temperature. A reliable lower limit can only be imposed by phenomenological arguments 


around the nucleosynthesis values, i.e. in the MeV range (Giudice et al. 2001). On the 
other hand, upper limits, taking into account the full spectrum of infiaton decay products 
in the thermalization process, are close to 10^^ GeV (McDonald, ^000 ). The main point 
for the antiproton emission study is that there is clearly a critical reheating tempera¬ 
ture, K. 10® GeV . If TfiH > the standard mass spectrum today is nearly not 

modified by the finite horizon size after inflation, whereas if T^h < the minimal 

mass becomes so large that the emission is strongly reduced. It is important to notice, 
first, that a rather high value of makes the existence of pbhs evaporating now into 
antiprotons quite unlikely, second, that the flux is varying very fast with T^h around 

rpc 

^RH' 

It could also be mentioned that above some critical temperature, the emitted Hawking 
radiation could interact with itself and form a nearly thermal photosphere (Heckler |l997|). 


This idea was numerically studied (Gline et al. 1999) with the full Boltzmann equation for 
the particle distribution and the Hawking law as a boundary condition at horizon. For the 
antiproton emission, the relevant process is the formation of a quark-gluon plasma which 
induces energy losses before hadronization. We have taken into account this possible 


effect using the spectrum given by Gline et al. ( 1999 ) : di^N/dEdt oc exp{—E/To) where 
Tq « 300 MeV normalized to the accurately computed Hawking spectrum. The resulting 
effect, dramatically reducing the evaluated flux, is shown in Fig. ||. 


3. Propagation model 

3.1. Features of the model 

The propagation of cosmic rays throughout the Galaxy is described with a two-zone 
effective diffusion model which has been thoroughly discussed elsewhere (Maurin et al. 
2001 - hereafter Paper I, Donato et al. 2001a - hereafter Paper H). We repeat here the 
main features of this diffusion model for the sake of completeness but we refer the reader 
to the above-mentioned papers for further details and justifications. 

The Milky-Way is pictured as a thin gaseous disc with radius i? = 20 kpc and 
thickness 2h = 200 pc (see figure where charged nuclei are accelerated and destroyed 
by collisions on the interstellar gas, yielding secondary cosmic rays. That thin ridge is 
sandwiched between two thick confinement layers of height L, called diffusion halo. The 
five parameters of this model are Kq, S, describing the diffusion coefficient K{E) = 
KoPTZ~^, the halo half-height T, the convective velocity 14 and the Alfven velocity 14. 
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Fig. 2. Primary antiproton flux with the standard mass spectrum before propagation (in 
arbitrary units) with different reheating temperatures. Curve (1) is for Trh = 3 • 10® 
GeV, curve (2) is for Trh = 10® GeV, curve (3) is for Trh = 6 • 10® GeV, curve (4) is 
for Trh = 3 • 10® GeV and the thick line is the spectrum without any cutoff. 


Specific treatment related to p interactions (elastic scattering, inelastic destruction) and 
more details can also be found in Paper II. 

Actually, a confident range for these five parameters has been obtained by the analysis 
of charged stable cosmic ray nuclei data (see Paper I). The selected parameters have been 
employed in Paper II to study the secondary antiproton flux, and are used again in this 
analysis (for specific considerations about the Alfven velocity, see section 5.2 of Paper II). 
In principle, this range could be further reduced using more precise data or considering 
different sorts of cosmic rays. For the particular case of /3 radioactive nuclei, Donato et 


al. 2001b showed that with existing data no definitive and strict conclusions can so far 


be drawn. We thus have chosen a conservative attitude and we do not apply any cut in 
our initial sets of parameters (which can be seen in figures 7 and 8 of Paper I). 
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Fig. 3. Primary antiproton flux with the standard mass spectrum. The higher points 
correspond to the pure Hawking spectrum and the lower points to the modified law with 
a QCD halo. 


3.2. Solution of the diffusion equation for anti protons 

Antiproton cosmic rays have been detected, and most of them were probably secondaries, 
i.e. they were produced by nuclear reactions of a proton or He cosmic ray (cr) nucleus 
impinging on interstellar (ism) hydrogen or helium atoms at rest. When energetic losses 
and gains are discarded, the secondary density satisfies the relation (see Paper H for 
details) 


2hS{z) q^^^{r,0,E) 


+ 


2hS{z)rf^ NP{r,0,E) 


NP{r,z,E) , 


(4) 


as long as steady state holds. Due to the cylindrical geometry of the problem, it is easier 
to extract solutions performing Bessel expansions of all quantities over the orthogonal 
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Fig. 4. Schematic view of the axi-symmetric diffusion model. Secondary antiproton 
sources originate from Cr/ism interaction in the disc only; primary sources are also dis¬ 
tributed in the dark halo which extends far beyond the diffusion halo. In the latter case, 
only sources embedded in the diffusion halo contribute to the signal (see Appendix P) . 


set of Bessel functions Jo{Qx) {Q stands for the ith zero of Jq and i = 1.. .oo). The 
solution of equation may be written as (see eqs A.3 and A.4 in paper II) 

= ^ 9 -(iJ)xexp||j^| |sinh|| (L-|.|)|/sinh||A|| , (5) 

where the quantities Si and Ai are defined as 

r t /2 a 2 'j 1/2 r S' I I 

5, = |^+4^| and A{E) = 2 + V, + K Smooth I . (6) 


We now turn to the primary production by pbhs. It is described by a source term 


distributed over all the dark matter halo (see section 2.1) - this has not to be confused 
with the diffusion halo - whose core has a typical size of a few kpc. At z = 0 where fluxes 
are measured, the corresponding density is given by (see Appendix ^) 


Arf'P"™(0) = exp 
where 


y^iL) 


-VcL 


2K J AiSinh(5*L/2) 




y.(L) = 2 / exp ( ^(L - z') ) sinh qr^{z')dz'. 




(7) 


( 8 ) 
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This is not the final word, as the antiproton spectrum is affected by energy losses when 
p interact with the galactic interstellar matter and by energy gains when reacceleration 
occurs. These energy changes are described by the integro-differential equation 

A,Nf + 2hdE {hLsi.E)Nf - Ki^{E)dENf} = 2h {qr"'iE)+qriE)+qriE)} (9) 

We added a source term ql‘^^{E), leading to the so-called tertiary component. It corre¬ 
sponds to inelastic but non-annihilating reactions of p on interstellar matter, as discussed 
in Paper II. The resolution of this equation proceeds as described in Appendix (A.2), 
(A.3) and Appendix (B) of Paper II, to which we refer for further details. The total 
antiproton flux is finally given by 

iVP'*°‘(i?o,0,A) = f2 + Jo(q^) (10) 

where IVf’*“(0, A) and Nf’^^^"^{0,E) are given by formulae (||), (^) and (^). We em¬ 
phasize that the code (and thus numerical procedures) used in this study is exactly the 
same as the one we used in our previous analysis (Paper I and II), with the new primary 
source term described above. 

As previously noticed, the dark halo extends far beyond the diffusion halo whereas its 
core is grossly embedded within L. We can wonder if the external sources not comprised 
in the diffusive halo significantly contribute to the amount of p reaching Earth. A careful 
analysis shows that in the situation studied here, this contribution can be safely neglected 
(see Appendix P). 


4. Top of atmosphere spectrum 

4.1. Summary of the inputs 

For the numerical results presented here, we have considered the following source term. 


which is a particular case of the profiles discussed in section 

El + Eq 


2.1 




m 


X (A) 


( 11 ) 


where (E) was shown in figures (0), (D and (^) for various assumptions. It corre¬ 

sponds to the isothermal case with a core radius Rc = i kpc (the results are modified 
only by a few percents if Rc is varied between 2 and 6 kpc). A possible flattening of this 
halo has been checked to be irrelevant. Actually the Moore or NEW halos do not need 
to be directly computed for this study as they only increase the total pbh density inside 
the diffusion volume for a given local density: as a result the flux is higher and the more 
conservative upper limit is given by the isothermal case. 

The primary source term (eq. 0) is injected in eqs. (0) and (0), then added to 
the standard secondary contribution (eq. 0), and propagated (eq. 0) for a given set of 
parameters (see section 3.1). Once this interstellar flux has been calculated, it must be 
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corrected for the effects of the solar wind, in order to be compared with the top of 
atmosphere observations. We have obtained the TOA fluxes following the usual force field 


approximation (see Donato et al. 2000 and references therein). Since we will compare our 
predictions with data taken during a period of minimal solar activity, we have fixed the 
solar modulation parameter to (j) = 500 MV. 

To summarize, the whole calculation is repeated for many input spectra, 

and all possible values of the propagation parameters 14, L, Kq, 14 S before being modu¬ 
lated. The modulated flux is then compared to data to give constraints on the local pbh 
density ■ 


4.2. Features of the propagated pbh spectrum 

Fig. I shows the top of atmosphere antiproton flux resulting from the standard mass spec¬ 
trum without any inflation cutoff {Thh > Trh) ^ local density = 10“^^ gcm“^. 

The two curves correspond to the extreme astrophysical models considered as acceptable 
in the extensive study of nuclei below Z = 30 (paper I). We can notice that uncertainties 
due to astrophysical parameters are very important on the primary propagated p flux. 
The degeneracy in the diffusion parameters that one can observe for stable nuclei - and 
also for secondary antiprotons (see figure 7 of Paper II) - is broken down for pbh p’s. 

This can be easily understood: secondary p and CR nuclei are created by spallations in 
the galactic disc, so that all the sets of diffusion parameters that make the primaries cross 
the same grammage (~ 20 g cm“^) give the same secondary flux. The simplest diffusion 


models (no reacceleration, no galactic wind ~ see for example Webber et al. 1992) predict 
that the only relevant parameter is then K^jL. On the other hand, primary antiprotons 
sources are located in the dark matter halo, and their flux is very sensitive to the total 
quantity of sources contained in the diffusion halo, i.e. to L. This explains the scatter 
of about one order of magnitude in predicted fluxes that are shown in Fig. for a given 
local density. 


4.3. Comparison with previous works 

In our models, all sets of parameters compatible with B/C data have been retained. In 
particular, a wide range of values for L is allowed (T > 1 kpc) and each L is correlated 
with the other parameters. The upper limit is set to 15 kpc, which is motivated by 
physical arguments (see for example Beuermann et al. 1985| and Han et al. 1997 by direct 


observation of radio halos in galaxies, showing that L is a few kpc, and Dogiel 1991 for 
a compilation of evidences). 

To our knowledge, the most complete studies about primary pbh antiprotons are 


those by Mitsui et al. ( 1996 ) and Maki et al. ( 199t ) in which the authors use roughly 
adjusted values for diffusion parameters and restrict their analysis to two cases, L = 2 
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PBH antiproton spectra 



Fig. 5. TOA primary antiproton flux with the standard mass spectrum after propagation 
for = 10“^^ gcm“^ and extreme astrophysical models. 


kpc and L = 4 kpc (and only one value of Kq/L = 0.008 kpc Myr“^). This is somewhat 
a crude estimation of what we know (or don’t know) about L. Moreover, in their models, 
neither galactic wind nor reacceleration is considered. The second point should not be 
very important (see Paper II) whereas the first point has a crucial impact on the detected 
antiproton flux. It can be shown that fluxes are at least exponentially decreased by the 
presence of a Galactic wind: setting 14 = 0 may overestimate the number of primary p 
reaching Earth. 

Effects of propagation parameters on primaries are qualitatively discussed in 


Bergstrom et al. ( 1999 ) for the case of SUSY primary antiprotons. At variance with all 
previous works on the subject, our treatment allows a systematic and quantitative esti¬ 
mation of these uncertainties, as the allowed range of propagation parameters is known 
from complementary cosmic ray analysis. 

The difference in the spectral properties we consider are less crucial: the input pbh 
spectrum used here is mostly the same as in Maki et al. ( 1996| ), whereas the accurate 
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dimensionless absorption probability for the emitted species is taken into account instead 
of its high energy limit and the possible presence of a QCD halo is also investigated. 


5. Upper limit on the pbh density 

5.1. Method 


Of course, the previously mentioned uncertainties will clearly weaken the usual upper 
limits on pbh density derived from antiprotons. To derive a reliable upper limit, and to 
account for asymmetric error bars in data, we define a generalized 


X 2-ii +a*'‘+ (Qi))^ oflP (tji) 3*, ) 


-E. 






+( 7 *' 


-(Qd)' 




where and o-“p+ (cr‘^“ and 17 “^“) are the theoretical and experimental positive 
(negative) uncertainties. Superimposed with BESS, CAPRICE and AMS data (Orito et 
al. 200C| , Maeno et al. 2000 , CAPRICE Coll. 2001 , Alcaraz et al. 2001 ), the full antiproton 
flux, including the secondary component and the primary component, is shown on Fig. ^ 
for 20 values of logarithmically spaced between 510“^^ and 10“^^ gcm“^. The 

standard mass spectrum is assumed and one astrophysical model is arbitrarily chosen, 
roughly corresponding to the average set of free parameters. 

An upper limit on the primary flux for each value of the magnetic halo thickness is 
computed. The theoretical errors included in the function come from nuclear physics 
{p + He^p+X and He + He ^ p+X, Paper II) and from the astrophysical parameters 
analysis which were added linearly, in order to remain conservative. The resulting x^ leads 
to very conservative results as it assumes that limits on the parameters correspond to 1 
sigma. 


5.2. Results 

Fig. ^gives the as a function of for L = 3 kpc. The horizontal lines correspond to 
63% and 99% confidence levels. In this paper the statistical significance of such numbers 
should be taken with care : they only refer to orders of magnitude. As expected, the 
value is constant for small pbh densities (only secondaries contribute to the flux) and is 
monotonically increasing without any minimum : this shows that no pbh (or any other 
primary source) term is needed to account for the observed antiproton flux. 

Fig. H gives the upper limits on the local density of pbhs as a function of L with 
the standard mass spectrum. It is a decreasing function of the halo size because a bigger 
diffusion region means a higher number of pbhs inside the magnetic zone for a given 
local density. Between L = 1 and L = 15 kpc (extreme astrophysical values), the 99% 
confidence level upper limit goes from 5.3 • 10“^^ gcm“^ to 5.1 • 10“^^ gcm“^. 
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PBH antiproton spectra 



Fig. 6. Experimental data from BESS95 (filled circles), BESS98 (circles), CAPRICE 
(triangles) and AMS (squares) superimposed with mean theoretical pbh spectra for 
between 5 • 10“^^ g.cm“^ (lower curve) and 10“^^ gcm“^ (upper curves). 


An important point shown in section (2.2.2) is that we are sensitive only to pbhs 
with masses between lO^^g and lO^'^g. The upper limit given here on the total local mass 
density of pbhs can therefore be given in a “safer” way as a numerical density of pbhs 
integrated between lO^^g and lO^^g. Although less interesting from the point of view 
of cosmology, this value has the great interest of being independent of the initial mass 
spectrum shape as this remains within the “increasing” part of the mass spectrum. The 
resulting numerical density is Uq pbh = 3.9 • 10“®^ cm“^ for L = 3 kpc. 

If we take into account the possible QCD photosphere (see figure ^) around pbhs, 
the previous upper limit is substantially weakened: < 1.2 • 10“^^ gcm“^ for the 

’’usual” L = 3 kpc. In this case, gamma-rays are a much more powerful tool to study 
PBHS. It should, nevertheless, be emphasized that this model is still controversial. First of 
all, because the quarks and gluons which are below threshold for pion production seem to 
be simply ignored, then, because the assumption that the Brehmsstrahlung interaction 
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Fig. 7. between the experimental data and the theoretical flux as a function of the 
local density of pbhs for a magnetic halo thickness of 3 kpc. The horizontal lines are 99% 
and 63% confidence levels. 

has a range of l/we could be wrong for hard interactions, which would lead to a great 
overestimation of the photosphere effect (Cline, private communication). 


6. Discussion and future prospects 

6.1. Comparison with other existing limits 


Many upper limits on the pbh explosion rate have been derived thanks to 100 MeV, 1 
TeV and 100 TeV 7 -rays or low-energy antiprotons. 

In the ultra-high energy range, a reliable search for short gamma-ray bursts radiation 
from an arbitrary direction have been performed using the CYGNUS air-shower array 
(Alexandreas et al. 1993| ). No strong one second burst was observed and the resulting 
upper limit is jdt d^V < 0.9-10® year“^pc“®. Very similar results were derived by 

the Tibet (Amenomori et al. 1995) and the AIROBIC collaborations (Funk et al. 1995). 
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Fig. 8. Upper limits on the local density of pbhs as a function of the magnetic halo 
thickness L. 


TeV gamma-rays have also been used to search for short time-scale coincidence events. 
The bursts detected are compatible with the expected background and the resulting 


upper limit obtained with 5 years of data (Connaughton 1998 ) is /dtd^V < 

3 • 10® year“^pc“®. 

The limit coming from antiprotons has been advocated to be far better: the previous 
study from Maki et al. ( 1996| ) gives dN^^J^/dtd^V < 2 • 10“^ year“^pc“®. However, we 
emphasize that this limit does not take into account the wide range of possible astrophys- 
ical uncertainties (in particular L, which can affect the limits by one order of magnitude). 
Moreover, we believe that the explosion rate is not the pertinent variable to use when 
comparing results from different approaches: thresholds differences between experiments 
make the meaning of ’’explosion” very different. With a mass spectrum oc for small 
masses, the number of exploding pbhs depends strongly on the value of the threshold. 
It makes the comparison of our results with the ones from Maki et al. very ambiguous. 
This is why we prefer to give our upper limit, either as a local mass density, assuming a 
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standard TUdss spectrum, < 5.3-10“^^ gcm“^, either as a number density, indepen¬ 
dant of the mass spectrum shape outside the relevant interval, < 1.3-10“^° cm~^ 

(whatever L). 

Gamma-rays in the 100 MeV region provide a sensitive probe to the presence of PBHs 
along the line of sight up to redshifts as large as ~ 700. Gamma-rays in this energy range 
have little interactions with the intergalactic medium and can travel on cosmological 
distances. The integration of the signal involves therefore a much larger scale than in the 
case of Milky Way antiprotons. It should also be pointed out that in this case most of the 
PBH population is involved as the dominant emission peaks above 100 MeV even at the 
present epoch. By matching the PBH cosmological emission to the extra-galactic gamma- 
ray diffuse background (MacGibbon & Carr 1991), the limit flpsH < 1-8 x 10~® h~‘^ was 
obtained. Being mostly based on the assumption that a standard PBH mass spectrum 
holds above M^,, this result is robust. Our limit {Pq^^ < 5.3 x 10“^^ gcm“^) also requires 
the same assumption. It no longer depends on the details of cosmic-ray propagation as 
it corresponds to the minimal possible value of 1 kpc for L. This result is therefore 
quite conservative. In order to discuss it in the light of the gamma-ray constraints, it 
should be noticed that any PBH population is a particular form of cold dark matter. 
When the latter collapses to form galactic halos and the intra-cluster medium, PBHs 
merely follow the collapse. Their abundance in the solar neighbourhood should trace their 
cosmological contribution to the overall value of ~ 0.3 — 0.4. Because a canonical 
isothermal halo has a solar density of ~ 0.3 GeV cm“^ - ~ 5.3 x 10“^® g cm“^ - we 
infer an upper limit of ~ 10“® on the contribution of PBHs to the galactic - and as 
mentionned above to the cosmological - dark matter. Our antiproton bound translates 
into flpBH < 10“® Hm ~ 4 X 10“®. Such a constraint is comparable to the limit derived 
from gamma-ray considerations. Antiprotons are produced by PBHs that are exploding at 
the present epoch and the limit which they provide is complementary to the gamma-ray 
constraint. 


6.2. Possible improvements on the antiproton limit 

New data on stable nuclei along with a better understanding of diffusion models (see 
for example Donato et al., in preparation) could allow refinement in propagation to 
constrain L. It is also important to notice that the colour emission treatment could 
probably be improved. Studying more accurately the colour field confinement and the 


effects of angular momentum quantization, it has been shown (Golubkov et al. 2000) 
that the meson emission is modified. Those ideas have not yet been applied to baryonic 
evaporation. 

Several improvements could be expected for the detection of antimatter from pbhs in 


the years to come. First, the AMS experiment (Barrau 2001) will allow, between 2005 and 
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2008, an extremely precise measurement of the antiproton spectrum. In the meanwhile, 


the BESS (Orito al. 2001) and PAMELA (Straulino et al. 2001) experiments should have 
gathered new high-quality data. The solar modulation effect should be taken into account 
more precisely to discriminate between primary and secondary antiprotons as the shape 
of each component will not be affected in the same way. The effect of polarity should 


also be included (Asaoka et al. 2001). 

To conclude, the antideuteron signal should also be studied as it could be the key- 
point to distinguish between PBH-induced and SUSY-induced antimatter in cosmic rays. 
Although the p emission due to the annihilation of supersymetric dark matter would 
have nearly the same spectral characteristics than the pbh evaporation signal, the an¬ 
tideuteron production should be very different as coalescence schemes usually considered 


(Chardonnet et al. 1997) cannot take place between successive pbh jets. 


Acknowledgements. The authors would like to thank J.H. Me Gibbon for very helpful discussions 
and D. Page for providing the absorption cross-sections. 


Appendix A: Computation of the primary galactic flux 

The solution for a generic source term z, E), is obtained by making the substi¬ 

tution 


2h6{z)q^^%r,0,E) < -> z, E) 

in equation i)- The Bessel expansion of this term now reads 

OO 

gP—(r,z,E) = J2qr"^iz,E)Jo(c.^) , 

i—1 

with (introducing p = r/R) 


(A.l) 


(A.2) 


(A.3) 


gf—(^, e;) = / p gP—(p,z, e;)Jo(c,p) dp 

The procedure to find the solutions of equation (^) is standard (see Paper I, Paper II 
and references therein). After Bessel expansion, the equation to be solved in the halo is 

dz^ Kdz i?2_ 

Using the boundary condition Nf{z = L) = 0, and assuming continuity between disc 
and halo, we obtain 

Ve{\z\-L)\ y,{L) 


K 


Aff’P"™(z) = -qf^'^iz) 


(A.4) 


Nf 


,prim / \ 

^ [z) = exp 


2K 


Ai sinh(5'iL/2) 
(Uc + 2 hr“^) 


cosh(S'iz/2) -I- ' ^ ^ sinh(S'i2;/2) 


Vijz) 

KS, 


where 


y,[z) = 2^ exp - z')^ sinh {^[z - 2 ')^ gf”’"(z')dz' 

In particular, at z = 0 where fluxes are measured, we have 

^ f -VcL \ _ y^{L) _ 

* ^ ^ 2K J AiSinh{S,L/2)' 


(A.5) 

(A.6) 

(A.7) 
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Appendix B: Sources located outside the diffusive halo 


By definition, diffusion is much less efficient outside of the diffusive halo, so that an¬ 
tiprotons coming from these outer sources propagate almost freely, until they reach the 
diffusion box boundary. At this point, they start to interact with the diffusive medium 
and after travelling a distance of the order of the mean free path, their propagation be¬ 
comes diffusive. Thus, the flux of incoming antiprotons gives rise to three thin sources 


located at the surfaces cti, (72 and as (see figure BT), the first with a ^ distribution rj^z) 
which has non-zero values only for z ^ L (the upper boundary of the diffusive volume), 
the second with a r distribution ri(r) which has non-zero values only for r R (the side 
boundary) and the third located at z « —L (the lower boundary). 


B.l. Solar density of antiprotons coming from external sources 

Let us focus first on the top and bottom surfaces ai and a^ located at z ~ L. As 
expression (j^) was obtained for a symmetrical source term, the corresponding external 
source term can be directly written as 

g“‘(r-, z) = r^‘^{r)p{L - z) with J p{z)dz = 1 

where r]{z)dz is the probability for an incoming nucleus to interact with the diffusive 
medium in the layer comprised between distances z and z -I- dz from the boundary, and 
jg antiproton flux coming from outside the diffusive halo. The function r]{z) is 
expressed in kpc“^. Using the fact that the surface source is located around small values 
of u = L — z, insertion of the above relation in formula (^) gives 

y^iL) « Y / ^V{u)du = J^r(r)S.u 

where is the Bessel transform of JF“'^(r) and where the mean free path u = 

J UT](u)du has been introduced, so that (^ becomes 

'-VcL\ 


Nr\0) = exp 


(B.l) 


2K J A,sinh(S',L/2)' 

Before turning to the computation of the incoming flux fF™°(r) at the box boundary, 
we can make a remark about the above expression. If galactic wind and spallations are 
neglected, it can be simplified into 


Af“‘(0) 


K cosh(^iL/i?) 

It occurs that the diffusion coefficient K and the mean free path u are related in a 
way that depends on the microscopic details of the diffusion process. For hard spheres 
diffusion, u = ZKjv so that 


Af“*(0) 


V cosh(((iL/i?) 
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The physical interpretation of this expression is interesting: for low Bessel indices i, i.e. 
for large spatial scales, the cosh(CiT/i?) term is close to unity and the Bessel terms in 
the disk density are proportional to those in the source term. For large Bessel indices i, 
i.e. for small spatial scales, the Bessel terms are exponentially lowered in the disk density 
compared to the source term. The small scale variations of the source are filtered out by 
the diffusion process. This filtering is more efficient for large values of the diffusion zone 
height L. 

In the following, we show the results for the full expression 


iV“‘(0) = 
where the 


(-VcL\ 

V 2K ) pcAi sinh(S',L/2) 
hard sphere expression has been used for u. 


(B.2) 



B.2. Computation of the incoming flux from external sources 

The flux of nuclei reaching the surface from a source located outside the diffusive halo is 
qpBH{r) dV 


= cos6» 




The integration over the angular coordinate (3 is performed using the property (see 

2E(a:) 


Gradshteyn 198C ) 
r d/3 


/o (a - 6cos/3)3/2 (a - b)^/a + b 


where E is the elliptic function (denoted Rp in Numerical Recipes, 1992) and where 
X = y/2bj{a + b). Finally, we can write 

= r r dPBHjr) _ 2E(a;) {z-L)RdRdz _ 

Jr=oJz=l 47r ((z - L)2 + {R- Roy) ^{z - Lf + (R + R^f ' 

This quantity is then numerically computed and Bessel-transformed to be incorporated 


in eq. (B.2). 
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B.3. Conclusion: the external contribution is negligible 


The contribution of these external sources is shown in table B.l for L = 1 kpc (the lowest 
value we allowed) and L = 5 kpc, with or without galactic wind. It is larger for lower 
halo heights L, but it is always less than 10“^ and it can be safely neglected. 


Table B.l. Fraction Next/Ntot of the antiproton density in the solar neighbourhood due 
to the external primary sources, for the pbh density profiles discussed in the text. The 
relevant cross-sections for antiproton spallations have been considered in the Ai term, 
and a kinetic energy of 1 GeV has been assumed. The results are not very sensitive to 
this particular energy value. 


density profile 

Vo = 

Okm/s 

Vo = 

lOkm/s 

L = 1 kpc 

L = 5 kpc 

L = 1 kpc 

L = 5 kpc 

Moore 

1.7 X lO""' 

4.9 X 10“® 

4.9 X 10"® 

4.0 X 10“^® 

NEW 

1.8 X lO""* 

5.4 X 10“® 

5.4 X 10"® 

4.3 X 10“^® 

Modified isothermal 

2.5 X lO""* 

9.3 X 10“® 

1.5 X 10"® 

7.0 X 10“^® 

Isothermal (a = 1 kpc) 

1.7 X 10"'^ 

4.3 X 10“® 

1.0 X 10"® 

3.6 X 10“^® 


Similar results apply to the side boundary <72 which is about Rgai — Rq ~ 12 kpc 
away. Roughly, the contribution is the same order of magnitude than would be obtained 
with a thin source cti located at z = 12 kpc, which can be neglected as discussed above. 

Note also that we present here only two combinations of 14 and L to derive the surface 
contribution, but all the good sets of propagation parameters used throughout this paper 
(and also in Paper I and Paper II) point toward the same values, which are not larger 
than those derived in table 1. We can thus conclude that primary surface contributions 
are negligible in all cases. 
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